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The semihadronic tau decay width allows a clean extraction of the strong coupling constant at 
low energies. We present a modification of the standard "contour improved" method based on a 
derivative expansion of the Adler function. The new approach has some advantages compared to 
contour improved perturbation theory. The renormalization scale dependence is weaker by more 
than a factor two and the last term of the expansion is reduced by about 10%, while the renormal- 
ization scheme dependence remains approximately equal. The extracted QCD coupling at the tau 
mass scale is by 2% lower than the "contour improved" value. We find as{M%) — 0.1211 ± 0.0010. 

PACS numbers: 13.35.Dx, 12.38.Bx 



I. INTRODUCTION 



Quantum Chromodynamics, QCD, is the theory of strong interactions, being the quark and gluon fields its basic 

degrees of freedom. It describes a rich variety of phenomena including confinement, chiral symmetry breaking, binding 

1—1, of hadrons and asymptotic freedom. For some observables a perturbative expansion in powers of QCD coupling Us 

i-G ' is possible, if the relevant energy scale is bigger than the QCD scale E'qcd ^ 1 GeV. The strong coupling is a 

'T^, fundamental parameter of the Standard Model and its determination is relevant in itself and for the identification of 

O ■ new physics. Extractions of the strong coupling constant from experiments covering energy scales from M^ = 1.78 

D \ GeV to ~ 200 GeV are consistent at a 1% level (at the Mz scale), providing an impressive test of asymptotic freedom 

r^ ^ • over an energy range of two orders of magnitude [1] . 

In this note we study the extraction of as from inclusive semihadronic tau decay or, more precisely, from the ratio 
^f) , of semihadronic to leptonic tau decay widths 

>: 

^: r. r(r^had;/, (7)) 

'NT I 

iy~j . This observable offers a unique and clean way of testing perturbative QCD (pQCD) at low energies, because of the 

f^ relatively large mass of the heaviest lepton that allows us to use pQCD and because of the inclusive character of the 

^D , observable avoiding the complication of hadronization. The theoretical expression for Rt is an integral in energy ^/s, 

^""^ ' from y^ = to -^/s = Mt, of the discontinuity of the two-point ly-channel correlation function, n(s), times a kernel. 

^ I Of course, due to the low energies involved and because n(s) is being evaluated at the physical cut, perturbative 

"'~j QCD cannot directly be used. However, thanks to the analytic properties of the exact n(s), Rr can be expressed 

rN as a s-plane contour integral along the circle \s\ = M^, with the Adler function in the integrand @t51. This allows 

» . a perturbative evaluation of n(s) and therefore of Rt, because the scale involved \s\ = M^ corresponds to a small 

" " ' absolute value of as and the contributions from the physical cut are suppressed (see details below). 

Experimentally we can separate final states with net strangeness from final states without net strangeness, and 
within the latter, the vector (V) from the axialvector (A) channels: 

Rr = Rf + R^ + R^. (2) 

If we are interested in the extraction of as from r decays it is convenient to use the vector plus axial part of 
^, RY^^ = RY + Rti instead of the V and A parts separately. There are two main reasons for this. Firstly, 
the experimental separation between V and A channels introduces an extra uncertainty not present in the sum, and 
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secondly, within the Operator Product Expansion (OPE) the experimentaUy extracted non-perturbative contribution, 
Spfp, and its associated uncertainty are suppressed in V+A compared with the V or A channels taken separately, due 
to a partial cancellation in the sum [a|- The experimental value of R^^^ is obtained from Rt and i?f . The completely 
inclusive hadronic quantity is obtained from the measured leptonic branching ratios 



Rr = ^ = ^ - 1-9726 == 3.640 ± 0.010, (3) 



and the updated strange part is 



Using these two values we have 



i?f = 0.1615 ±0.0040. (4) 



i?^+^ = 3.479 ±0.011. (5) 

This is the main experimental input for the extraction of the strong coupling constant. The extracted value of the 
MS coupling constant from R^^^ in Contour-Improved perturbation theory (CIPT) at N'^LO order is 

af{Ml) = 0.1217 ± 0.0017, (6) 

after renormalization group (RG) evolution up to Zo scale. We observe a tension when comparing with the world 
average obtained by Bethke [l] 

Rcthkc/ 11 #2 



a. 



(M|) =0.1184 ±0.0007, (7) 



dominated by the lattice QCD extraction from the HPQCD collaboration [3| 



, lattice/ 71 r2 



a= 



(M|) =0.1183 ±0.0008. (8) 



On the other hand, when comparing with the value extracted from Zq decays, also at N'^LO order [g|, 

af (M|) = 0.1190 ± 0.0026, (9) 

or with the value extracted from the pt dependence of the inclusive jet cross section in pp collisions at ^/s — 1.96 
TeV, the most precise determination from hadron-hadron colliders [9|, 

af{Ml) = 0.116llO;S°48 , (10) 

we observe agreement within the uncertainties (with a higher alpha central value from tau decays). 

After the publication Q of the af correction io Rj, several groups have used this result for the extraction of a^ 
and condensates in the context of the OPE [g, [^ Il0l - [l3j . The values of as obtained from these groups are partly 
incompatible with each other [1|, mainly due to the usage of either CIPT or Fixed Order Perturbation Theory (FOPT). 
The groups also have differences in the way they treat non-perturbative contributions. 

Alternatively, evaluation of low energy QCD observables can be performed in the context of analytic QCD [ij] 
(anQCD; for reviews and further references see [15^-17] ). Such approaches have a number of advantages. For example, 
the running coupling af'XQ'^) is an analytic function of Q^ [Q^ = —q^), without the Landau singularities outside 
the negative semiaxis in the complex Q^-plane. Therefore, evaluated expressions for n(s = — Q^) in anQCD lead to 
identical results for the integrals eqns. ([TT|) and p^ . This is not true in perturbative QCD, due to Landau singularities 
of as (Q^) on the positive Q^-axis. Some anQCD models [la [3 can reproduce the correct values of i?^+^ by adjusting 
some low-energy free parameters, but then some other attractive features are lost. However, other anQCD models, 
with otherwise attractive features, have a tendency to give too low values of i?^+"^, cf. Refs. [ISi [2(|422| . 

The main point of this work is to present an alternative evaluation of Rt within pQCD. We call this procedure 
modified CIPT. Due to the relatively small energy involved, there are important effects coming from the manner we 



use the renormalization group, as we know from the difference between CIPT and FOPT. It turns out that in modified 
CIPT the extracted value of as{M^) is 2% lower than in CIPT, decreasing the difference between the value obtained 
from tan decays and the world average from [l|. Using as input the experimental values RY~^^ (eq. ([5])) and Snp (see 
below), we obtain in modified CIPT a'^'^^uf) = 0.341 ± 0.008, and in CIPT ^^(M^) = 0.347± 0.015. Evolving the 
coupling up to the Zq scale we get, respectively, af'^^M^) = 0.1211 ± 0.0010 and aP(M|) = 0.1217 ± 0.0017. The 
quoted uncertainties are total uncertainties. 

The modification of CIPT is simple. Instead of the usual power series expansion a + ci a^ + C2 a^ + . . . (where: 
a = tts/Tr), the Adler function is expressed by a nonpower series of the form a + ci 02 + £2 03 + . . ., where dn+i are 
proportional to the n'th derivative of the coupling a{Q'^) and c„ are the new expansion coefficients. Thus, the terms 
of the series are proportional to the coupling and its derivatives, i.e. proportional to the coupling, the /3-function and 
derivatives of the /3-function. Therefore, it can be said that the /3-function plays in modified CIPT a more central 
role than in CIPT. This expansion in derivatives of as was introduced in [19| in the context of skeleton-motivated 
expansion and analytic QCD. It turns out that compared to CIPT the new expansion shows a lower RG dependence 
and hence a lower theoretical error within the method. 

The total hadronic ratio in e^ e~ collisions, i?e+e-('S), is another timelike observable which can be expressed in 
terms of the corresponding Adler function. Using the new expansion for the Adler function we present a new and 
simpler expression for the RG-improved i?e+e-(s) in terms of the new couplings a„. 

In Section In] we review the standard evaluation of Rt- Section IlIII contains the main part of this note, here we 
present and study the new approach. The uncertainty in the extraction of a^ from i?^+ is discussed in Section IIV[ 
while in Section |V] a new expansion for the RG-improved Re+e- (s) is presented. Finally, the conclusions are given in 
Section EH 



II. CIPT AND FOPT 

The semihadronic tau decay ratio can be expressed as |23] 



''^^ ds / s \V. . s \1 



^^ = io m['-m) ['^'mln''^''^'^^ ^''^ 



where n(s) is the correlator of two T^-channel currents (for more detailed expressions see the reviews [2J, [2 

n(.) = K,\'{nU^) + ii^A^)) + Ks\'{iiUs) + iitis)), (12) 

with Vij the elements of the Cabibbo-Kobayashi-Maskawa matrix. Eq. ([TT|) cannot be evaluated directly using 
perturbative QCD due to the small energy involved in the integral. However there is a way out. By general arguments 
we know that the exact function n(s) is analytic function in the whole s-complex plane, excluding the physical region 
s > 0. Therefore, using the Cauchy theorem we have 



|s|=Af2 

and integrating by parts 

1 /■ (7t 1 

i?,. = — i -il^a:)'{l + a:)-D{~xM^), (14) 

zm J X 2 

\x\ = l 

where the Adler function D{Q^), defined by 

is an observable (the renormalization scheme dependent constant of n(s) is eliminated). Note that in eqns. ([T5)) and 
(ITil) the quantities n(s) and D{Q^) are evaluated at an absolute value of energy square equal to M^. Therefore the 



absolute value of the complex running coupling constant is small enough in order to perform a perturbative treatment. 
The theoretical expression for the vector plus axial part of Rr can be written in the Operator Product Expansion 
framework as 

RX+^ = 3|Kd|'5e^(l + 5^ + C + '^2 + 5np), (16) 

where the perturbative QCD correction, the central object of this article, is given by Jq. Sew = 1.0198 ± 0.0006 |26l | 
and 5'ew = 0.001 ± 0.001 [27| are electroweak corrections, 62 = (—4.3 ± 2.0) x 10^^ are light quark masses effects, 
5jvp = (-5.9 ± 1.4) X 10-3[i] are non-perturbative contributions,^ and Kd = 0.97418 ± 0.00027 [28ij. From these 
values and eq. ([5]) we obtain the experimental value 

(5o = 0.204 ± 0.004, (17) 

the number from which we extract as using only pQCD. 
The perturbative QCD contribution can be expressed as 

(5o = ^ y —{l~xf{l+x)D{~xMl), (18) 

\x\ = l 

where the reduced (canonically normalized) Adler function D{Q^), defined from eqns. ([T?l) . (ITS)) , and p^ . is the 
massless QCD perturbative (leading twist) contribution to 



3|KdP5, 



ew 



Perturbatively and without considering quark masses the vector and axialvector contributions are equal. The vector 
(or axial) Adler function is known at fourth order in QCD 

4 n-l 

D{Q') - Y. «"(^') E ^»." log"(gVM'), (19) 

where a{Q^) = a{Q'^)/Tr, /i is the renormalization scale, and Cn^m the expansion coefficients. Only the coefficients 
Cnfl are independent. By using the renormalization-group (RG) equation, the coefficients Cn^m with to > 1 can be 
obtained as linear combinations of the Cn',o with n' < n, with coefficients depending on the perturbative /?- function. 
The Adler function is itself an observable and its RG improved expression is 



^''''(Q') = E^».o«"W')- (20) 

The /3-function coefficients are normalized as 

a, ° , = p{a) = -(/3o a" + A a^ + p^ «* ^ ^^ ^5)_ (21) 

a log ^^ 

In Contour Improved Perturbation Theory (CIPT) the expression for (5o, eq. ([TS]). is evaluated using expression ([20)) 
for the Adler function, where the behavior of the coupling a(Q^) is exactly dictated by eq. (|^T)) . Thus, the Adler 
function is RG improved along the contour (|T5)) . which is the right thing to do considering that d{Q^), and therefore 
the integrand of eq. ()18p . is an observable. An alternative procedure is Fixed Order Perturbation Theory (FOPT) 



^ This value is obtained in CIPT. We assume that <5jvp does not vary in modified CIPT. 
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FIG. 1: Left: Real part of a2(M^e' ) (solid line) as a function of 9 compared to o?(M^e'^^) (dashed line). In both cases we take 
a(M^) = 0.340/7r. Right: The corresponding imaginary part. 



where eq. ([TO)) is used in (fT5|) . choosing a unique renormahzation scale /i^ — M^ along the contour, 
values of the coupling extracted using these two methods deviate significantly from each other 



The central 



„ci 



a 



FO 



(M^) 
(M^) 



0.347, 
0.326. 



(22) 
(23) 



Another argument in favor of CIPT is the fact that this procedure is much more stable under renormahzation scale 
variations than FOPT Q.^ We will not consider FOPT in the following. From the discrepancy between the two 
methods we learn that in the extraction of a^ from i?T- there is an important effect coming from the manner the RG 
is used. This motivates the next section. 

Before presenting the derivative expansion, let us make a further comment. It is well known that near the 
Minkowskian semiaxis the perturbative Adler function (or n(s)) does a bad job at reproducing the exact func- 
tion due to the resonance structure. Therefore the factor (1 — x)^ in eq. ([TS)) is crucial for a clean evaluation of Rr, 
because it suppresses the contribution from this problematic region.^ 



III. MODIFIED CIPT 



In this, the central part of the article, a new method for the evaluation of Rr is presented. We call the new approach 
modified CIPT. As in CIPT, the semihadronic tau decay ratio is evaluated in eq. (IT51) using for the Adler function 
D a RG-improved expression. The proposed modification consists in using instead of the standard series in powers of 
a{Q'^) given in eq. (l20l) a nonpower expansion for the Adler function in terms of the new couplings dn{Q^)- Truncated 
at the last known term the new expansion is given by 



^(Q') = I]c„a„(Q2), 



(24) 



with the tilde couplings defined as 



^ For discussions and comparisons between both approaches see [a. llOl. [29ll3C|. 
^ Duality violations are studied in l3ll. |32| . 



0-m+l 



(-1)" d"'a 
13^ m\ d(logQ2)r 



(25) 



The derivatives are evaluated perturbatively using eq. ((2T|) . The new coefficients c„ are obtained from the coefficients 
Cm,o with m < n. The tilde couplings are normalized such that dn = a" + ©(a"^^). Note that in the new expansion 
5,1 = a, 0,2 = — /3(a)/^o, ^3 = /^('3')/3'(a)/(2/3o)i 6tc. The beta function and its derivatives play in ([M|) a more central 
role than in eq. (j20p . The real and imaginary parts of the ffist three new couplings, 02, 0,3 and 04, are plotted along 
the integration circle of eq. (|18p in Figs. [1] [2]and[3l respectively, together with the corresponding powers a". Note 
that in all three cases dn{M^) > a"{M^). In general, the ratio dn/a" growths with n: 



^ = l + r„ + 0(a2), 

1 "' 



(26) 



where 



'n+l 



r„(n + l)/n + A/^o- 



(27) 



Only for low values of a and n we have dn/a" near 1. The analyticity properties of the tilde and standard couplings 
are similar. If we restrict ourselves to real Q^ , there are poles and cuts in the infrared region. 

The series (UHl) and 1^^ are formally equal if infinite number of terms in both expressions is considered. However, 
it is believed that these series are asymptotic. If we had the complete perturbative series we would truncate them at 
their respective (in absolute value) smallest terms in order to give a meaning to the sums. In each case, the difference 
between the sum and the exact value of the original quantity is expected to be smaller than the last considered term. 
As a consequence, the rearrangement or reshuffiing of eq. (|20l) in eq. (l24t is not immaterial: it leads in general to 
different values of the truncated series for the Adler function. Furthermore, with the presently known coefficients the 
asymptotic behavior of the Adler function is possibly not reached and we are forced to truncate the series before, 
obtaining also different values in GIFT and modified GIFT. 

Next we evaluate Rr and compare GIFT and modified GIFT. The MS three-flavor Adler function coefficients c„_o 
and c„ are 



ci,o = 1, C2,o = 1.6398, c3,o = 6.3710, C4,o = 49.076, 
5i = 1, £2 = 1.6398, 53 = 3.4558, £4 = 26.385. 



(28) 
(29) 



The first two coefficients are equal by construction and the third and fourth tilde coefficients are about half the 
standard ones. The value of the first unknown coefficient has been estimated [a] using the method of "fastest apparent 
convergence" to be cs^o = 275, corresponding to C5 = —25.4.* 
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Z2n = l 


ELi 


a 




0.1132 
0.1082 


0.0210 
0.0244 


0.0092 
0.0081 


0.0081 
0.0107 


(0.0051) 
(-0.0019) 


0.1515 
0.1515 


(0.1566) 
(0.1496) 


0.3556/7r 
0.3400/7r 



TABLE I: Contributions to the Adler function using the standard power expansion (first row) and the tilde expansion (second 
row). In both cases the coupling a is chosen in order to obtain D{M^) = 0.1515 when summing up to n = 4. Numbers in 
brackets consider the estimation for the first unknown coefScient, 05,0 = 275. 



As an intermediate step, we first study the Adler function at the scale Mr using the power series eq. (f20| and the 
tilde expansion eq. ([M)) . The results are shown in Table H] We demand in both cases D{M^) — 0.1515 when summing 
up to n = 4, extracting a different value of strong coupling a(Af^) in each case. We observe in modified GIFT a value 
of the coupling 4 — 5% smaller than in GIFT, this difference is much smaller than the experimental error for the Adler 
function 33]. The small ratio c„/c„ g for n = 3 and 4 is compensated at this scale by a high ratio o„/a". If we take 
the estimate 05,0 — 275 (see the numbers in brackets) then this term is suppressed in the tilde expansion. 



** If the unknown coefficient is estimated as cs^o = C4,o(c4,o/c3,o) = 378, we obtain £5 = -1-77.6 
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FIG. 2: Left: Real part of a3{M^e'^ ) (solid line) as a function of 9 compared to a^(M^e'^^) (dashed line). In both cases we take 
a(M^) = 0.340/7r. Right: The corresponding imaginary part. 

Note also that the n = 4 contribution is close to the n = 3 one. This fact could be accidental or can be interpreted 
as the term n = 3 (or 4) being the smallest term of the asymptotic series for the Adler function. In the latter case, 
for the best estimation of the Adler function we must truncate the series at this term. However, the tilde expansion 
for the Adler function is thought as an intermediate step in the evaluation of Rt or related quantities. 
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Z^n = l 


EL-, 
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CI 

c/ 


0.1513 
0.1484 


0.0308 
0.0372 


0.0128 
0.0104 


0.0090 
0.0078 


(0.0038) 
(-0.0001) 


0.2038 
0.2039 


(0.2077) 
(0.2037) 


0.347/7r 
0.341/7r 



TABLE II: Contributions to the semihadronic tau decay width using CIPT (first row) and modified CIPT (second row). In 
both cases the coupling a{M-r) is chosen in order to obtain So = 0.204 when summing up to n = 4. Numbers in brackets 
consider the estimation for the first unknown coefficient, C5,o = 275. 

We evaluate 5q in MS scheme using the expression (Ull), eq. ^^ for CIPT (CI), and eq. (HI) for modified CIPT 
(C/). The results are shown in Table [TTl Again, the value of a{Mr) is taken different in the two approaches in order 
to obtain the same value of Rr, when summing up to n = 4, Sq — 0.204. In modified CIPT we get a value of the strong 
coupling at the r mass scale by about 2% lower than in CIPT. In addition, we observe in modified CIPT compared 
to CIPT a smaller last term of the series (n = 4 term). We see no signal of having reached the smallest term of the 
asymptotic series for R^. Including the estimate C5.0 = 275 the modified CIPT series has a surprisingly low last term. 
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So, CI 
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0.3831/7r 
0.3400/7r 
0.2812/7r 


0.2009 
0.1984 
0.1907 


0.2020 
0.2031 
0.1991 



TABLE III: Renormalization scale dependence of Sq. The absolute value of the (squared) renormalization scale is chosen to be 
/i^ = ^Mr with 5 = 0.7, 1, and 2. We take as reference a{M^) = 0.34/7r. 



We study the renormalization scale dependence of So in MS scheme, comparing CIPT and modified CIPT. The 
(squared) renormalization scale is chosen to be fi^ — ^M^ with ^ = 0.7, 1, and 2, and we use as reference a{M^) — 
0.34/7r. The results for ^o are shown in Table IIIII Taking as a measure of the scale dependence of 60 its range of 
variation when ^ varies between 0.7 and 2, we obtain 0.0102 and 0.0040, for CIPT and modified CIPT, respectively. 
If we extract a^ from the experimental value of 5q, these uncertainties translate to Aa{M^) =0.013 and 0.005 for 
CIPT and modified CIPT, respectively. Thus, using this criterion, renormalization scale dependence is by more than 
factor two weaker in modified CIPT than in the standard CIPT. 
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FIG. 3: Left: Real part of a4{M^e ) (solid line) as a function of 6 compared to a^(M'^e*'^) (dashed line). In both cases we take 
a{M'^) = 0.340/7r. Right: The corresponding imaginary part. 
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TABLE IV: Renomialization scheme dependence of So. The quantity So is evaluated in the 't Hooft scheme extracting a(tH) 
in CIPT and modified CIPT. In the last column the value of a(tH) in converted to the MS scheme. 



Finally, we study the renormalization scheme dependence in the extraction of a^ from (5o, comparing MS and 't 
Hooft (tH) schemes. The latter scheme is defined by setting all but the first two coefficients of the beta function 
equal zero, i.e. (32 — Ps = ■ ■ ■ = 0. The result of the extraction of a in the 't Hooft scheme is shown in Table IIVI 
Comparing these values with the values of (|ll|, we obtain Aa(M^) = 0.004 in both cases. Thus, renormalization 
scheme dependence is approximately equal in CIPT and modified CIPT (it is slightly weaker in modified CIPT) . 

Modified CIPT is a new kind of perturbative expansion. The approach is valid in itself and possesses some attractive 
properties as is a lower renormalization scheme dependence. However, from the point of view of the standard power 
series of the Adler function, mCIPT performs the sum ([20| up to n = 8 considering nonzero coefficients c„,o for n = 5 
to 8. For example, when expanding the truncated expression (j24p in powers of a, we obtain cs^o = 300.4. The natural 
question to ask is how good this value compares to the exact one. To answer it we need to calculate the corresponding 
Feynman diagrams, a task not expected to be done in the near future. What we can do as a test of the method is to 
compare its prediction for the known coefficients C3_o and 04^0- From ci_o and C2,o the estimate for 03^0 is 2.92, and 
from ci^Oi C2.0 and C3.0 the estimate for C4.0 is 22.7. Comparing with the exact values, cf. eq. (j28p . we see that in 
both cases mCIPT includes a significant part of the next term of the power series. Therefore, at least in these two 
cases, mCIPT is an improvement also from this point of view. Both estimates (for c^^ and 04^0) are lower than the 
exact value by a factor 2.2 (2.19 and 2.16 respectively). 



IV. UNCERTAINTY IN THE EXTRACTION OF a. 



The experimental uncertainty in the extraction of as, Aa'^'^P = ±0.005, comes from the uncertainty in the value of 
the pseudo-observable quantity 60 given in eq. (|17l) . By far the main contribution here is due to the experimental 
value of R}^~^"^, given in eq. ([5]). The non-perturbative contribution Snp is considered here as an input for So and 
also contributes to the experimental uncertainty of ag. Its value and its associated uncertainty, which are obtained 
from the moments of Rt, are rather small [6j] and the uncertainty in Snp has almost no relevance for the uncertainty 
of So (while the effect of Snp on the central value of So is 0.006, i.e. 1.5 times So's experimental uncertainty). 

The theoretical uncertainty, within modified CIPT, was obtained in the previous section varying the renormalization 
scale, Aa'''^' = 0.005, and scheme, Aa^'^^ = 0.004. Adding them in quadrature we obtain Aa^^°° = 0.006. Then, the 



extracted value of the strong coupling constant at the M^ scale in modified CIPT is 



af^^iMl) = 0.341 ± 0.005'='^P ± 0.006"^™, (30) 

= 0.341 ±0.008. 

For comparison we give the corresponding values in CIPT: 

af{M^) = 0.347 ±0.005<=''P±0.014*'''=°, (31) 

== 0.347 ±0.015. 

Alternatively, we could estimate the theoretical uncertainty for a^ as coming uniquely from the way we use the 
RG: taken the difference between the extracted central values using CIPT and modified CIPT. Coincidentally, this 
would lead to the same theoretical uncertainty obtained above within modified CIPT. We are not allowed to sum 
them. This would imply a double counting because, by definition, the difference between CI and modified CI is given 
by higher order contributions. 

Conventionally, we compare the values of the strong coupling extracted from different experiments at a particular 
scale, the Mz scale. We perform the evolution at four loops, with three-loop matching conditions [3J| at the thresholds 
/ithr — '2,mq (q — c,b). The uncertainty ±0.0005 due to RGB evolution is as given in 6]. Evolving eq. (15(1)) from Mt 
to AIz we get 

a^^^(M|) = 0.1211 ±0.0006'=''P±0.0007*'''=°±0.0005'=™\ (32) 

= 0.1211 ±0.0010, 

in modified CIPT. 

V. HADRONIC RATIO IN e+ e" COLLISIONS 

Another Adler function related quantity is the ratio of the total hadronic to muonic cross sections, i?g+e-(s), 
which is proportional to Imn(s) (with n(s) in the neutral channel). Analogous to the definition of D{Q^) a reduced 
Re+e- (s), fi(s), can be defined. The function R{s) can be written as a function of the reduced Adler function as 

— s+ie 

The non-RG-improved expressions for R{s) and D{Q^) differ in the so-called 7r^-terms, which are numerically impor- 
tant. A RG-improved expression for R{s) cannot be obtained directly from the non-RG-improved R{s), because the 
RG is valid in the Euclidean region and s is timelike. The right way to proceed is to apply the RG group to the Adler 
function, obtaining D^'^{Q'^), and then get an improved version of R{s) from eq. ([55]) [35]. 

Using the tilde expansion we obtain a simple expression for R{s), because the integral (j33p can be performed for 
all but the first term of the tilde series of the Adler function: 

-s+ie ^ 

m^T^ I -a(z) + ^-^^^ImK_i(-.-ze)}. (34) 

— s — ie 

Thus, we obtain a new expression for R^+e- (s)- The first term of eq. (|M)) is the Minkowskian coupling [3a], and higher 
terms of the series are proportional to the imaginary part of the tilde couplings evaluated at the time momentum s, 
times a suppressed coefficient. It would be interesting to perform a phenomenological analysis of Re+e- (*) using this 
new expression. 
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VI. CONCLUSIONS 

In this note we present a modification of Contour-Improved Perturbation Tlieory, i.e., of the standard method for 
the evaluation of the semihadronic tau decay width, R^, within the context of pQCD and the OPE. Due to the low 
energy involved the truncated af result is sensible to higher order terms. The way we use the renormalization group in 
calculating Rt from the Adler function leads to important uncertainties in the evaluation of Rr , as is well known from 
the difference between CIPT and FOPT. Both, in CIPT and in the proposed approach, the Adler function is evaluated 
using a varying renormalization scale and not a fixed one fi — Mr as in FOPT. The new ingredient in modified CIPT 
is in the series we use to evaluate the Adler function: instead of the usual power series expansion a + ci a^ + C2 a"^ + . . . 
(where a = as/ir), the Adler function is expressed by a nonpower series of the form a + ci 0,2 + C2 03 + . . ., where the 
new tilde couplings dn+i are proportional to the n'th derivative of the coupling a{Q^) and c„ are the new expansion 
coefficients. It can be said that the /3-function plays in modified CIPT a more central role than in CIPT. This 
expansion in derivatives of as was introduced in |19| in the context of skeleton-motivated expansion and analytic 
QCD. 

Modified CIPT has some advantages compared to contour improved perturbation theory. The renormalization scale 
dependence is weaker by more than a factor two and the last term of the expansion is reduced by about 10%, while 
the renormalization scheme dependence remains approximately equal. 

The total hadronic ratio in e'^ e~ collisions, i?e+e-(s), is another timelike observable which can be expressed in 
terms of the corresponding Adler function. Using the new expansion for the Adler function we present a new and 
simpler expression for the RG-improved i?e+e- (s) in terms of the new couplings a„. 

The extracted value of ag from the vector plus axial non-strange R^ is in modified CIPT 1.8% (0.5%) lower than 
in CIPT, at the r (Zq) scale. We obtain a'f^\M^) = 0.341 ± 0.008 and af'^^(M|) = 0.1211 ± 0.0010. 
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